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We study and compare equilibrium and aging dynamics on both sides of the ideal glass transition 
temperature Tmct- In the context of a mean field model, we observe that all dynamical behaviors 
are determined by the energy distance e to threshold - i.e. marginally stable - states. We furthermore 
show the striking result that after eliminating age and temperature at the benefit of e, the scaling 
behaviors above and below Tmct are identical, reconciling en passant the mean field results with 
experimental observations. In the vicinity of the transition, we show that there is an exact mapping 
between equilibrium dynamics and aging dynamics. This leads to very natural interpretations and 
quantitative predictions for several remarkable features of aging dynamics : waiting time-temperature 
superposition, interrupted aging, dynamical heterogeneity. 



Introduction 

Glassy systems are characterized by very slow relaxa- 
tions and non-equilibrium behavior. In general, there is 
however a high temperature (or low density, pressure, 
etc) equilibrium phase where all quantities relax to their 
equilibrium value in finite time. Between the two is the 
putative glass transition, which has been identified as a 
true singularity in the free energy in spin glasses or in the 
dynamic free energy in mean field models of glass formers. 
The main predictions of mean field models concerning the 
glass phase are the existence of an effective temperature 
Tf.f f at which the slowest degrees of freedom are "ther- 
malized" and a specific form for the decay of the cor- 
relators at long times. Although effective temperatures 
have been reported in experiments and numerical simu- 
lations [3] , the time evolution predicted by mean field 
models is in general not observed. There is thus still a 
deep challenge in interpreting the measurements of slow 
relaxations in the glass phase as well as in describing 
the crossover between aging and equilibrium dynamics 
which take place when one is patient enough to wait for 
the system to reach equilibrium. Among many interesting 
features of glassy materials, waiting time-temperature su- 
perposition (WTTS) has been reported in various glass 
formers |11|5J|5]. 

The aim of this Letter is, in the context of a mean field 
model, to answer the following important questions. Are 
mean field predictions compatible with experimental ob- 
servations ? Are slow dynamics in the glass phase intrin- 
sically different from slow dynamics in the equilibrium 
phase ? How does the crossover from aging to equilibrium 
occur at the critical point ? What is the origin of the ob- 
served WTTS ? The answer to these questions we will 
shed new light on the central role played by the relaxa- 
tion of the energy. 

The Letter is organized as follows : we first describe the 
model and the old solution of its dynamics, and then im- 
prove it with some new results. Next, we establish a cor- 
respondence between slow dynamics at equilibrium close 



to Tmct and aging dynamics below Tmct, which be- 
comes an exact mapping at Tmct- The consequences of 
this correspondence, WTTS, interrupted aging and dy- 
namical heterogeneity, are then discussed. Some results 
presented here follow from long and technical derivations, 
which will not be given here and will be detailed somew- 
here else 0. 



THE MEAN FIELD PICTURE OF AGING 

We start discussing the dynamics of the spherical p- 
spin model, where exact statements can be made [T]. 
It can be seen as a ground for establishing rigorous 
results within the landscape approach following Gold- 
stein [8], and developed by many authors - see [9l [10] 
for instance. The model consists of N variables Si 
subject to the constraint Sf — N, with Hamil- 
tonian H — V- - Ji, ... i Si, ■ ■ ■ Si , where the J's 
are quenched Gaussian random variables. Explicit equa- 
tions can be written for the dynamics of correlators 
C{t,t') = {Si{t)Si{t')) and response functions R{t,t') = 
{SSi{t) /Shi{t')) to an external field. The solution of these 
equations, which can be found in [1 , captures the fea- 
tures of glassy dynamics, where two well separated times 
scale emerge, corresponding to different kinds of degrees 
of freedom, some - "local" - responding as in equilibrium, 
and others - "structural" - being out of equilibrium. In 
the equilibrium time sector, fluctuation dissipation theo- 
rem (FDT) holds. In the long time aging sector, FDT 
does not hold, but instead a modified form of it, where 
the bath temperature T has to be replaced by an effective 
temperature T^f f > T which can be computed explicitly. 
A fine analysis shows that slow degrees of freedom are in- 
deed thermalized at T^f / rather than T, coming from the 
contribution of the structural degrees of freedom to the 
entropy at some energy Exh- The threshold energy Exh 
is the energy separating unstable and stable states, and 
corresponds to the asymptotic energy : at long times, a 
glassy system is close to marginal stability. In the equili- 
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brium phase, the energy is Eeq{T) is above ExhiTMCT), 
and at the approach of Tmct, the system gets closer to 
the threshold and slows down, almost trapped by mar- 
ginally stable states. Similarly, in the glass phase, where 
the equilibrium energy is below the threshold, the sys- 
tem loses equilibrium when it gets close to ETh{T) and 
thus ages. In both cases, the closer the system is from 
the threshold level, the slower it is [11]. Before making 
this more precise, we briefly recall the known scalings for 
the correlator C{t, t^) and the response function R(t, t^), 
using as usual t = t — , the equilibrium and glass 
phases [Ij. 

The old results 

In the equilibrium phase, i.e. T — Tmct + there 
are two characteristic time scales, the /3-decay around 
the plateau and the a-relaxation [T^] : — e~^/^° 

and Ta{e) — e^'', 7 = l/2a -I- 1/2&, where a and b are 
well known MCT exponents. There, the scalings for C 
are, in order : C{t, t^) = qc + ^c^ij jTjj) and C(i, t^,) = 
1cCa{T/Ta), whcrc is the non-ergodicity parameter at 
Tmct- 

In the glass phase T < Tmct, two characteristic time 
sectors have also been identified [T]. The first time sec- 
tor is stationary : when t ^ i^, there is time trans- 
lation invariance (TTI) : C{t,t^) = Cst(t) and FDT 
is still vaKd : R{t,tw) = The second - aging - 

time sector corresponds to r ~ t^, [TB]. There, FDT must 
be modified : R{t,t^) — 9C(t,t^) ^j^^ scahng is 
C(t,tjj^) — C{h{tyj)/h{t)). An interesting but unpleasant 
feature is that the analysis of the aging equations at infi- 
nite tw is not sufficient to determine h{t) [T], which has re- 
mained unknown until recently. From this point of view, 
equilibrium and aging dynamics seem to differ strongly. 
We now improve an analysis made recently of the dyna- 
mical equations at large but finite tw |16| and show that 
in fact both phases are very similar. 

New results about the aging regime 

Recently, it was shown that for T < Tmct, aging ac- 
tually sets in around the plateau [TSJ [T7] , a time sector 
which is usually just considered as the matching point 
of the TTI and aging sectors. It was shown that in this 
regime, the correlators have a scaling form C{t,t^) = 
Q + tw^QiT/^w), where the exponents are related through 
a = (3a. In addition, matching with the aging regime, 
the uncertainty on the function h(t) was considerably 

reduced to h{t) = exp (^Aj^^, where A is a constant 
which can be safely absorbed into the relaxation time 
and set to 1, while the aging exponent fi < I verifies 
fib = a + bp and now a and b are aging MCT expo- 



nents. Remark that /i < 1 is required for this solution 
to be consistent. Situations where this occurs have been 
reported in various experiments, e.g. [18]. Values of /i ex- 
ceeding 1 (super-aging, see e.g. [TH]) are the signature 
of a different relaxation mechanism. After this analysis, 
the complete solution is still to be found, as the expo- 
nents are unknown although all determined by a. Ho- 
wever, interestingly, in the same time the energy was 
shown to relax as E{t) = ETh+E2t~'^°' , meaning that the 
whole slowness of the correlators is actually encoded in 
the energy relaxation towards the threshold level. This is 
remarkable and will lead to the striking correspondence 
between equilibrium and aging dynamics, which we shall 
establish below. Last, but not least, the scaling of the res- 
ponse function is R{t,t^) = -t-"-l'^w'{T/t'^)/T and the 
ratio X{y) — w{y)/g{y) interpolates smoothly between 
X{-oo) = 1 and X{oo) ^ X ^ T/T^ff, showing that all 
interesting features of the aging regime already occur in 
the plateau regime, which corresponds to the time scales 
where motion propagates from local to structural degrees 
of freedom. 

AGING IS LIKE EQUILIBRIUM 

We now analyze further and improve these results, in 
order to obtain a unified picture of equilibrium and glass 
phases. First, we recall the simplest of the two equations 
verified by the scaling functions g and w around the pla- 
teau [TB] (where Xq is such that g(a;o) = 0) : 

^(x)+rdy^^^=^(.o), 

^p{x) = w{x)'^ + / dyw'{y){w{x — y) — w{x)) 
Jo 

When T = Tmct, T = Teff, and Eq. ^ reduces to the 
equation for the scaling function of equilibrium dyna- 
mics [12]. Second, writing h{t) — exp(j){t), it is easy to 
verify that in order to have C(h{ty^) /h{t)) be of order 1, 
one needs that r ~ (j)'{tw) = H^^- This has an important 
consequence : the scaling C{h{tw) /h{t)) with a stretched 
exponential h(t) is in general very hard to observe ex- 
cept in the p-spin model [17 , and in the large limit 
it differs only through tiny asymptotic corrections to the 
observed scaling C{t,tw) = C{T/tl^)- This explains the 
many reported discrepancies between mean field predic- 
tions and measurements [21 [3]. Third, and more strikin- 
gly, the aging relaxation can also be related to the scaling 
function of the a-regime in the equilibrium phase, at all 
temperatures. Indeed, the scaling function C was found 
in [I4 to verify 

(p - 1)(1 - ,)C(A) - C(\)'-'(l -1 + iX). 
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We now show that the solution of this equation is in 
fact a familiar function of the equilibrium phase. Writing 
A = exp(— r), one gets, after a bit of algebra : 

m + ^f^ f dtvp-\t)v'{T-t) = Q, (3) 

'^^MCT -If) 

with 'D{x) — C(exp(— a;)), which is exactly the equation 
verified by Cq [TH[T3]. Using h{t) « exp(T/tJ^), and re- 
marking that the small r behavior must match the (3- 
regime, we can determine unambigously T), leading to : 

C{t,t^) = qCa {T/Ta{t^)), whcrC Ta{t 



ages. This ends when the energy reaches its equilibrium 
value E{T), which gives : t^ross \T-Tmct\~^^^°'^ ■ This 
leads to the following scaling form for the a-relaxation 
time : t„ = T^^^{T)T{t^/tcross), where T{x) ~ a;'''^ at 
X 1 and T(oo) — 1. Remark that tcross may also be 
expressed as |r — Tmct\'^^^'' , that is the exponent is the 
ratio of the ones of the relaxation time in both phases, 
which gives an easier way to measure it in the MCT re- 
gime. 

Dynamical heterogeneity 



Mapping equilibrium onto aging 

The interpretation of the preceding paragraph will be 
central result. We introduce e(r, t^) — E{t^) — ExhiT) 
and obtain for T < Tmct ■ rpit^) = e(T,t^)-'^/'^°' and 
T„(t„) = e(T, t„)-i/2a-i/2b^ Added to the above scahngs, 
this leads us to the conclusion that once temperature and 
dependences are expressed in terms of e{T,tw) only, 
equilibrium and aging dynamics almost map onto each 
other. Almost, because a, b and the scaling functions in 
the /3-regime explicitly depend on T below Tmct- Ho- 
wever, the mapping becomes exact at the transition, and 
is very accurate close to it (see Fig. [T]). This is a rather 
unexpected and spectacular result, as it gives a route 
to studying aging dynamics from the knowledge of the 
deeply supercooled equilibrium phase - calculations in 
both phases are not intrinsically different. In addition, 
it incorporates in a natural way the WTTS reported in 
experiments. 



Interrupted aging 

An immediate consequence of that is interrupted aging. 
Indeed, the glass transition temperature Tg is where a 
glass former seems to lose equilibrium, but where, with 
a bit of patience, it is still possible to have it equih- 
brate. More precisely, there is a crossover waiting time 
tcross such that for t^ < tcross aging occurs, while for 
tw > tcross, the system has reached equilibrium and its 
dynamics is TTI. By definition of Tg, tcross is large. A 
naive guess for the value of tcross would be tcross = , 
the equilibrium relaxation time. However, this is too 
naive, as this is the time for the system to have visited 
several equilibrium states, while tcross is the time where 
aging stops and is in principe much smaller. In our model, 
it is possible to compute tcross using the previous analy- 
sis, from the following gedanken experiment : quench the 
system from T = 2Tmct down to T = Tmct + £, e < 
During the first stage of the dynamics, the energy is close 
but not enough to its equilibrium value and thus relaxes 
as a power law : E{t^) — ETh{Td) « t~'^°''' and the system 



Using the previous results, it is possible to describe 
the general behavior of correlation functions generally 
used for describing dynamical heterogeneities (DH). 
We start with the now widely studied X4(^7^io) = 
/ dy (^{SCix, t, t^)6Cix + y, t, t^)) - {5C{x, t, t^)f^ , 

where 5C{x,t,tw) is a local correlator at position x. It is 
possible [7] to extend the field-theoretic approach to Xi 
used by Biroli and Bouchaud [20 to the aging regime. It 
results that XA{t,t^) is a sum of ladder diagrams built 
with the full correlator 5C{x,t,tw), which time depen- 
dence follows that of the previous paragraphs. This has 
several immediate consequences. First, all scaling results 
given in previous analysis of Xi ™- the equilibrium phase 
can be directly applied to the glass phase, provided 
we replace e by e(T, f^), providing MCT predictions 
for DH in he glassy phase. Second, it was reported 
in [21 that, once X4(^7 iii;)/niaxx4(i, f^) is plotted versus 

r 

1 — C{t,t^), all data - i.e. for all t^'s - collapse onto 
a single master curve. This could be obtained from 
general scaling considerations, but interestingly, it comes 
out naturally from the fact that all tu,-dependance of 
both quantities plotted is through e{T,t^), which gives 
the natural parameterization of the master curve. In 
addition, close to Tmct, the full master curve should 
coincide with the one of the equilibrium phase. Another 
quantity which has been shown to be interesting above 
Tmct is Xt{t) = ^'^Jf ^ , which has been shown to have 
the same critical behavior as X4(t), and which is more 
easily accessed experimentally 22J. It is very simple to 
generalize this quantity to the glass phase, as d/dT is 
identical to d/dE, up to the non-singular multiplying 
factor OE/dT. It is thus natural to generalize xt to : 



}Cw ('^ tw : tw ) 



dC{T + tw,tw) 

de{T,tw) 



(4) 



In the equilibrium phase, Xw{t,tw) reduces to 
Xt{t), while in the glass phase, it becomes 

( ^'^Ot'*"'"' ) '^'"^"^ot*""*™'' ■ T^^^s expression makes Xw 
as easy to measure in experiments as xt- In Fig. [T] Xw 
at T = Tmct + 10"^ and xt at T = .5 w Tmct - 0.125, 
with approximately the same value of e are plotted 
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depend on temperature, and thus the predictions made 
here should in general be tested at close temperatures. 
This would also be of valuable help for determining si- 
tuations where the energy landscape approach gives a 
qualitatively correct picture. On this prospect, and more 
speculatively, one may ask whether the effect of activa- 
tion at low temperature may be also reabsorbed in the 
same way. It would be also worth coming back to data 
from previous experiments, where disagreement with for- 
mer mean field predictions where found, and reinterpret 
them using the predictions given in this Letter. 



Fig. 1: Comparison of xt{t) above Tmct (straight line) and 
Xw{t) well below Tmct (dashed line), with the same value 
of e. The Curves have been rescaled vertically, as the factor 
de{T,tw) / dtw occurring in Xw is not known. 

versus t. Clearly, the peaks of the a-relaxation coincide 
very well, while the /3-regime slightly differ, having 
two different exponents, respectively a « 0.396 and 
a w 0.448. When the same procedure is made at 
T = Tmct for Xw, the curves fall on top of each other. 
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CONCLUSION 

In this Letter, we have shown that, within a mean 
field model, dynamics in the ideal glass phase are essen- 
tially identical to equilibrium dynamics just above the 
glass transition, and that all important features in both 
phases may be absorbed in a single function e{T,tw), 
with important physical significance, being the distance 
in energy to threshold states. This is a remarkable result, 
as it opens the door to interpreting the aging data ob- 
tained in the glass phase of molecular or colloidal glass 
formers. In particular, it provides general scaling laws for 
multipoint correlation functions - relaxation, dynamical 
heterogeneities - as well as the natural parameterization 
to seek, e(T, tw), which may alternatively be an enthalpy 
difference. It also naturally predicts WTTS, which was 
reported in aging measurements On the conceptual 
aspect, it is very unexpected, because it shows that aging 
dynamics is actually very similar from equilibrium dy- 
namics in the deeply supercooled regime, which is ra- 
ther far from common thinking about aging. But this is 
not so surprising, and reflects the fact that both in li- 
quid and glass phases, slow dynamics occurs because of 
the roughness of the energy landscape, which shape does 
not change qualitatively more than static quantities when 
crossing the glass transition. Here, several strong predic- 
tions have been made, which should be tested in expe- 
riments and numerical simulations. Doing so, one must 
keep in mind that in general, the exponents a and h both 
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